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Abstract 

The enumeration of lattice paths in wedges poses unique mathematical challenges. 
These models are not translationally invariant, and the absence of this symmetry com- 
plicates both the derivation of a functional recurrence for the generating function, and 
solving for it. In this paper we consider a model of partially directed walks from the origin 
in the square lattice confined to both a symmetric wedge defined by Y = ±pX, and an 
asymmetric wedge defined by the lines Y = pX and Y = 0, where p > is an integer. We 
prove that the growth constant for all these models is equal to 1 + v2, independent of the 
angle of the wedge. We derive functional recursions for both models, and obtain explicit 
expressions for the generating functions when p = 1. From these we find asymptotic 
formulas for the number of partially directed paths of length n in a wedge when p = 1. 

The functional recurrences are solved by a variation of the kernel method, which we 
call the "iterated kernel method". This method appears to be similar to the obstinate 
kernel method used by Bousquet-Melou (see, for example, references [5, 6]). This method 
requires us to consider iterated compositions of the roots of the kernel. These compositions 
turn out to be surprisingly tractable, and we are able to find simple explicit expressions 
for them. However, in spite of this, the generating functions turn out to be similar in 
form to Jacobi ^-functions, and have natural boundaries on the unit circle. 

PACS numbers: 05.50.+q, 02.10.Ab, 05.40.Fb, 82.35.-x 
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1 Introduction 



The problem of counting random walks on lattices with various restrictions is perhaps one of 
the oldest problems in enumerative combinatorics, with a history that dates back at least 100 
years [1]. It has also seen a great deal of recent activity, particularly surrounding problems of 
counting random walks on the slit-plane and quarter-plane [15, 22, 13, 12, 4, 7, 5]. 

These models pose interesting mathematical problems, and powerful methods have been 
developed in recent years to solve for the generating functions of path problems. Normally, 
these methods are a three step process: First a recurrence is determined, this is solved in the 
second step, and lastly, the asymptotics of the number of paths are extracted. 

Perhaps the simplest and most studied model is Dyck paths. While there are numerous 
techniques for enumerating Dyck paths, the most powerful technique involves a recurrence for 
the generating function which is solved and expanded to determine an explicit expression. If 
d n is the number of Dyck paths of half-length n, then the generating function g t = ^ n>0 d n t n 
satisfies the recurrence 



so that d n is given by Catalan's number. 

In this paper we follow a similar strategy to determine the generating function and asymp- 
totic expressions for the number of partially directed paths confined to a wedge. The wedge 
destroys translational invariance in the model, and both the derivation of a recurrence for the 
generating function, and solving the generating function, poses difficult mathematical problems. 
Our strategy is in principle no different from the above for Dyck paths - we shall derive func- 
tional recurrences for the generating functions, solve those in special cases, and then find the 
asymptotics for the number of paths. Unfortunately, the problem for general wedges appears 
intractable, and even in the cases that we do solve we encountered significant difficulties. 

Models of paths and walks frequently appear as simple models of polymers in dilute solution 
in the physics literature [29]. The properties of polymers are in part determined by their 
conformational entropy, and models of walks and paths contributes to our understanding of the 
significance of the conformational entropy contributions in the free energy of polymers. These 
entropic contributions are important when polymers are in confined geometries. For example, 
the steric stabilisation of colloids by polymers results when polymers are confined to the spaces 
between colloidal particles [21]. This situation have been modelled by studying paths confined 
to the slab between two planes, see for example [9, 23]. 

Lattice random walk models of polymers in confined geometries are generally more tractable. 
These models can generally be solved, at least in principle, by a Bethe ansatz or constant term 
formulation. This technique have been used to solve for random walks in a half-space and 
which interacts with the boundary of the space [20]. Such random walk models, however, do 
not take into account the volume exclusion of monomers in a polymer. A more realistic model 
is the self-avoiding walk [19]. This model is non-Markovian, and while much is known about it 
from constructive [17, 18] and conformal invariance techniques (in two dimensions) [8], solving 
it remains beyond the current techniques in combinatorics. 

Self-avoiding walk models of polymers in confined spaces have not been solved (except in the 
most trivial of cases), but there are some results in the literature. For example, the exponential 
growth constant of self-avoiding walks in a wedge geometry is independent of the angle of the 
wedge [16]. Additionally, for self-avoiding walks in wedges, conformal field theory have been 
used to examine the dependence of scaling exponents on the wedge angle [8, 11]. 



gt = i + tg 2 t 



(i.i) 



with solution 
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Figure 1: Models of directed walks above a line Y = +pX. The walks are constrained to take 
only north and east steps (left); north, east and north-east steps (middle); and north, east and 
south steps (right). 

The introduction of directedness in a self-avoiding walk in a wedge may give models of 
directed or partially directed walks in a wedge which are both self-avoiding and which may in 
some cases be solvable. Models of directed and partially directed walks in wedge geometries (see 
Figure 1) have been studied previously [10, 28, 25]. These models include directed paths in a 
wedge; a model which is related to Dyck paths. It is interesting that the radius of convergence of 
the generating function is known in this model, even for wedges with wedge-angles of irrational 
cotangent [28, 25]. 

In this paper we consider models of a partially directed path confined in wedges (see Fig- 
ure 2). These models are similar to the directed path models in Figure 1, however, they are 
also substantially more challenging, since the path interacts with the wedge on two sides, rather 
than on only one side. As a result, it is much harder to find their generating functions and 
analyse their asymptotics. 

1.1 Directed and partially directed paths 

A directed walk on the square lattice is a path taking unit steps only in the north and east 
directions. Such objects are necessarily self-avoiding; they cannot revisit the same vertex. 
Partially directed paths may take unit steps only in the north, south and east directions with 
the further condition that no vertex is visited twice — ie they are self-avoiding. Hence, north 
steps cannot be followed by south steps and vice-versa. The generating function of such walks 
can be derived using standard techniques: 

W(t) = ^c^= (1.3) 

n>0 

where c n is the number of walks of length n and t is the length generating variable. An expansion 
of W{t) of W(t) in t produces an explicit expression for c n : 

c n = l - ((1 + V2) n+1 + (1 - v/2) n+1 ) • (1.4) 

The exponential growth constant is the exponential rate at which c n increases with n. This is 
given by 

fi = lim 4 /n = 1 + V2. (1.5) 

n— »oo 

This is the most fundamental quantity in this model from a statistical mechanics point of view. 
The radius of convergence of W(t) is /i" 1 , and the limiting free energy is the logarithm of 



3 




Figure 2: (left) The symmetric model of a partially directed path in a p-wedge formed by the 
lines Y = ±pX. (right) The asymmetric model of a partially directed path in a p-wedge formed 
by the line Y = +pX and the X-axis. 



the growth constant, k = log/x, this defines the explicit connection between the combinatorial 
properties of the model and its thermodynamic properties. 

We show some models of directed and partially directed paths in wedge geometries in 
Figure 1. The model in Figure 1 (left) was considered in [10, 27]. In general the growth 
constant is a (non-trivial) function of the wedge angle. The derivative of the free-energy with 
respect to the wedge angle gives the moment of the entropic force exerted by the polymer on 
the wedge and this was computed in [27]. This model may be also generalised by introducing 
an interaction between the line Y = +pX and the path, or by considering partially directed 
paths or Motzkin paths instead [26, 25]. 

In Figure 1 (right) a partially directed path confined to the wedge above the line Y = pX 
and the y-axis is proposed instead. This model was considered in reference [26]. If the partially 
directed path is instead confined to the wedge between the X-axis and the line Y = pX, then 
the model in Figure 2 (right) is obtained, which is the subject of this paper. 

In particular, we consider the variants illustrated in Figure 2 - firstly a model of a partially 
directed path in a wedge formed by the lines Y = ±pX (we call this the symmetric model - see 
Figure 2 (left)), and secondly a model of a partially directed path in a wedge formed by the 
X-axis and the line Y = +pX (this is the asymmetric model - see Figure 2 (right)). 

The related model of a partially directed path in a wedge with last vertex in the line y = pX 
is illustrated in Figure 3. This is a bargraph path above the line Y = +pX. This model was 
examined in reference [26], and while the generating function g p (t) is not known explicitly, it 
is given by 

*w = T^whm (L6) 

where h(t) is an appropriate solution of the equation 

Mt ) = ^ Ml + Mf)f ( 1 + T _m_), (1 . 7) 

where (as above), t is conjugate to the number of edges in the path. 
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Figure 3: A bargraph path above the line Y = +pX. The generating function of this model 
is not known explicitly. However, a set of equations derived in reference [26] can be solved to 
determine the radius of convergence of the generating function for integer values of p. 



The model in Figure 3 was also used as a model of adsorbing bargraphs which interact with 
the line Y = +pX [10], and an asymptotic expression for the adsorption critical point has been 
estimated in reference [26] (the location of the singular point on the radius of convergence of 
the generating function). 



1.2 Partially Directed Paths in Wedges 

Consider the square lattice Z 2 of points in the plane with integer coordinates. Let p > be an 
integer. The symmetric p-wedge V p is defined by 

V p = {(n, m) E Z 2 | where n > and —pn < m < pn} . (1-8) 

The (asymmetric) p-wedge is defined by 

W p = {(n, m) £ Z 2 | where n > and < m < pn} . (1-9) 

Let v Pt n (resp. w p ^ n ) be the number of partially directed walks in V p (resp. W p ) of length n. 

In the next section we establish some basic facts about the asymptotic growth of v n>p and 
w ntP as n — > oo. In Section 3 we find functional equations satisfied by the corresponding 
generating functions, which we solve in Sections 4 and 5. We then analyse the generating 
functions to determine v n ,\ and w n) \ to leading order. We show in particular that 

v n>1 = 4,(1 + V2) n + ,/ = (4 + (-l) n 4 + 0(l/n)) (1.10) 

for the number of paths in a symmetric wedge when p — 1 where A , 4 an d 4 are constants. 
The asymmetric wedge poses more difficult mathematical problems, and we were only able to 
show that 

:i + V2) n 



w n>1 = y r J-4 (Bo + o(l)) , (1.11) 
Vn + l 



for some constant Bq. 



2 Growth constants in wedges 

In this section we prove that the growth constant for partially directed walks is independent of 
the angle of the wedge and is equal to that of unrestricted partially directed walks. First, let b n 
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be the number of partially directed walks in the wedge defined by the lines X = and Y — 0, 
whose last vertex lies in the line Y — 0. These paths are counted by the generating function 
go(t) defined above, and singularity analysis gives the following lemma. 

Lemma 2.1. The growth constant of partially directed paths in the wedge defined by X = 
and Y = is 

lim b\i n = (1 + y/2) = At. (2.1) 

n— >oo 

This result can be used to determine the growth constants of partially directed walks in the 
wedges V p and W p . We first prove existence of the growth constants. 

Lemma 2.2. For any given p G (0, oo) the following limits exist: 

lim^ = ^ and lim w 1 ^ = fi. (2.2) 

The limits satisfy 

1% <^<^=(1 + V r (2)) (2.3) 

Proof. We have that w UtP < v UjP < c n . Hence, if the above limits exist, we must have aC — 

fi;<fi = (i + V2). 

To show existence, we prove that the sequences are super-multiplicative. Take any walk 
counted by v n>p and append a horizontal step, and any walk counted by v m ^ p . This gives a walk 
of n + m + 1 steps that lies within V p , and so is counted by v n+m+ i )P . Hence v ntP v mjP < v n+m+ i tP . 
A standard result (Fekete's lemma) on super-additive sequences (which we can apply by taking 
logarithms) then implies that ju£ exists. The proof for walks in W p is identical. □ 

Next, we show that — fj, p , and we show that they are equal to 1 + v2. 

Lemma 2.3. For any given p G (0, oo) we have 

b n < Wtf np -\ +nN+N ) )P (2.4) 

And hence lim n ^ 00 b\[ n < [i™. 

Proof. Take any walk counted by b n . By prepending \np] + 1 horizontal steps, this walk will fit 
inside the wedge, W p . Now append another horizontal step and a walk counted by b n — repeat 
this until there are iV walks counted by b n . This gives a walk counted by w^ np -\ +n N+N), P - Thus 
we have the first inequality. Taking logs and dividing by ( \np] + nN + N) gives 

N 1 

(\np] + nN + N) l ° gK ~ (\np]+nN + N) l0gW ^ +nN+N ^ ^ 

Take the limit as N — > oo to obtain 

-log6 n <log^. (2.6) 

Next, take the limit as n — > oo to complete the proof. □ 

By combining the above lemmas we can prove that the growth constant for partially directed 
paths is independent of the wedge angle. 

Theorem 2.4. For any given p G (0, oo) 

/ i; = / i; = /i = (i + v / 2). (2.7) 

This shows that the dominant asymptotic behaviour of the number of walks is independent of 
the wedge angle. Below, we show that the leading sub-dominant behaviour is also independent 
of the wedge angle [ie for p > 1). 
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3 Functional equations for walks in wedges 



3.1 The symmetric wedge model 

Consider a model of partially directed paths in a symmetric wedge as illustrated in Fig- 
ure 2 (right). If p is an integer or a rational number, then the path may touch vertices in 
the lines Y = ±pX. These vertices are visits in the lines Y = ±pX. In the event that p is an 
irrational number such visits cannot occur, however the path may approach arbitrarily close 
to the adsorbing lines (for large enough AT-ordinate) . In this paper we shall only consider the 
simplest version of this model, and we assume that p is a positive integer. Even in this case the 
model is apparently intractable, and we have only found the generating functions when p = 1. 

We will derive a functional equation satisfied by the generating function of partially directed 
paths in V p (those illustrated in Figure 2 (left)), by finding a recursive construction, similar to 
those in [3, 5] (and elsewhere). 

Let x be the generating variable for horizontal edges in the path and let y be the generating 
variable for vertical edges in the path. Introduce generating variables a and b to be conjugate to 
the distances between the last vertex in the path and the line Y = —pX and the line Y = +pX 
respectively. The generating function of the paths are now denoted by g p (a,b; x,y) = g p (a,b) 
where the variables x and y are suppressed. 

It turns out that the construction and resulting functional equation is simplified by consid- 
ering only those partially directed walks that are either a single vertex (no edges) or end in 
a horizontal step. Let f p (a,b;x,y) = f p (a,b) be the generating function of such paths. It is 
simply related to g p (a, b) via 

f p (a,b) = l + x(ab) p g p (a,b). (3.1) 

We now obtain a functional equation satisfied by f p by recursively constructing the paths 
column-by-column. Each path is either a single vertex, or can be constructed from a shorter 
path by appending either a horizontal step, or a sequence of up steps followed by a horizontal 
step, or a sequence of down steps followed by a horizontal step. See Figure 4. 
Consider a path counted by f p (a, b), and see Figure 4. 

• Appending a single horizontal step to its end increases the distance of the end point from 
both wedge boundary lines by p. Hence the generating function of paths with a horizontal 
edge appended is x(ab) p f p (a,b). 

• Appending an up step to the end of such a path increases the number of vertical steps by 
1, increases the distance from the line Y = —pX by 1 and decreases the distances from 
the line Y = +pX by 1. Hence such a path has generating function y(b/a)f p (a, b). Hence 
appending some positive number of up steps gives 1 ^° a /p(a, b). Appending a horizontal 

step to the end of such a path gives (by the above reasoning) x(ab) p 1 yb ^ a f p (a, b). 

• Similarly appending some positive number of down steps followed by a horizontal step 
gives x(ab) p ^j^f p (a, b). 

Unfortunately, when appending up or down steps it is possible that the resulting path will 
step outside of the wedge. Hence we must subtract off the contributions from such paths 
(Figure 4 right-top and -bottom). 

• Consider a path that ends at a distance h + from the line Y = +pX. It we append more 
than h + up steps to the path then it will leave the wedge. We can decompose the resulting 
path into the original path with exactly h + up steps appended, and an "overhanging" V 
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Figure 4: Constructing partially directed walks in the wedge V p . Every walk is either a single 
vertex, or can be obtained from a shorter walk by appending a horizontal edge (left), or a run 
of north steps and a horizontal edge or a run of south steps and a horizontal edge (centre-top 
and -bottom). Care must be taken to not step outside the wedge when appending north or 
south steps (right-top and -bottom). 



shaped path which is a sequence of some positive number of up steps and a horizontal 
step (see Figure 4 top-right). 

Appending exactly h + up steps to the path increases the distance from Y = —pX by 
h + , decreases the distance from Y = —pX to zero. This gives the generating function 
f p (by, b). The overhanging piece is (by the reasoning above) enumerated by x(ab) p 1 ^° a - 

Hence the g.f. of walks that leave the wedge is given by x(ab) p 1 _ f(py, b). 

• Similarly when appending too many down steps we obtain configurations counted by 

x ( ab ) p T^7bfp( a ' a y)- 

Using the above construction we arrive at the following theorem 
Proposition 3.1. The generating function f p (a,b;x,y) = f p (a,b) of partially directed walks 
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ending in a horizontal step in the wedge V p satisfies the following functional equation: 



f p (a,b) = l + x(ab) p f p (a,b) 




(f p (a,b)-f p (by,b)) 



(f P (a,b) - f P (a,ay)) 



(3.2) 



The generating function of all partially directed walks in V p is given by 



g p (a,b)=x 1 (ab) p (f p (a,b) - 1) 



(3.3) 



In the next section we turn to the problem of solving this functional equation. 
3.2 The asymmetric wedge model 

Let us now turn to the construction of partially directed paths in the asymmetric wedge W p . Let 
the generating function of all partially directed walks in this wedge be denoted k p (a, b; x, y) = 
k p (a, b) where the variables x and y are suppressed. 

As above, the resulting functional equation satisfied by the generating function is simpler 
if we consider only those walks that are either a single vertex or end in a horizontal step. Let 
this generating function be denoted h p (a, b; x, y). This is simply related back to k p by 



We now use the same construction as was used above for the symmetric case — each walk 
is either a single vertex, or can be constructed from a shorter walk by appending either a 
horizontal step, or a run of up steps and a horizontal step, or a run of down steps and a 
horizontal step — see Figure 5. Again care must be taken not to step outside the wedge, and 
so those walks that do step outside the wedge must be removed. Indeed the argument is almost 
identical to that used above, except that a horizontal step contributes xa p instead of x(ab) p , 
since a horizontal step increases the distance from the line Y = +pX by p, but does not change 
the distance from the line Y = 0. 

The above construction gives the following theorem: 

Proposition 3.2. The generating function h p (a,b; x,y) = f p (a,b) of partially directed walks 
ending in a horizontal step in the wedge W p satisfies the following functional equation: 




(3.4) 



h p (a, b) = 1 + xa p h p (a, b) 




(3.5) 



The generating function of all partially directed walks in V p is given by 




(3.6) 



We solve this equation in Section 5. 
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Figure 5: Constructing partially directed walks in the asymmetric wedge W p . Every walk is 
either a single vertex, or can be obtained from a shorter walk by appending a horizontal edge 
(left), or a run of north steps and a horizontal edge or a run of south steps and a horizontal edge 
(centre-top and -bottom). Care must be taken to not step outside the wedge when appending 
north or south steps (right-top and -bottom). 

4 Solving the symmetric case 

At first sight, one might try to solve equation (3.2) by the iteration method used in [3], however 
the coefficients of the equation are singular when a = by and b = ay. Multiplying both sides of 
the equation by (a — by) (b — ay) gives a non-singular equation, however when we set a = by or 
b = ay the equation reduces to a tautology. 

Instead we apply a variation of the kernel method, which we call the iterated kernel method. 
This appears to be similar in flavour to the "obstinate kernel method" used by Bousquet-Melou 
[5, 6]. We start by collecting all the f(a, b) terms together on the left-hand side of the equation 
— this gives the kernel form of the equation: 

K(a, b) f p (a, b) = X(a, b) + Y(a, b) f p (a, ya) + Z(a, b) f p (yb, b), (4.1) 
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where the functions K(a, b), X(a, b), Y(a, b) and Z(a, b) are given by 

K(a, b) = (b- ya)(a - yb)(l - x(ab) p ) - xy(ab) p (a 2 + b 2 - 2yab), (4.2a) 

X(a, b) = (b- ya) (a - yb) (4.2b) 

Y(a, b) = -xya p+1 b p (a - yb) (4.2c) 

Z(a, b) = -xya p b p+1 (b - ya) (4.2d) 

for each integer p > 1. The function K(a, b) is called the kernel of the equation. Note that the 
equation is symmetric under interchange of a and b: 

f p (a,b) = f p (b,a) K(a,b) = K(b,a) X(a,b) = X(b,a) Y(a,b) = Z(b,a). (4.3) 

We solve equation (4.1) by substituting an infinite number of pairs of a and b values that set 
the kernel K(a, b) to zero. While the method we describe below should work for general p, the 
resulting expressions are so complex that the process becomes intractable. This will become 
clear even in the case that p — 1, which we simplified only after significant effort. 



4.1 Iterated kernel method for V\ 

When p = 1, the kernel becomes a quadratic function of a and b and we can explicitly write 
down (simple) expressions for its zeros. This is not generally true for larger values of p, and 
since simplifying our expressions requires that compositions of the zeros of the kernel must 
simplify as well - the general p case appears intractable from a practical point of view. 

Thus, we restrict ourselves to p = 1. We write j\(a,b) = F(a,b) and the coefficients in 
equation (4.1) become 

K(a, b) = (xy 2 a 2 - xa 2 - y)b 2 + (1 + y 2 )ab - ya 2 (4.4a) 
X(a, b) = (b- ya)(a - yb) (4.4b) 
Y(a, b) = -xya 2 b(a - yb) (4.4c) 
Z(a, b) = -xyab 2 (b - ya) (4 Ad) 

Let /3±(a; x, y) = (3±(a) be the zeros of K(a, b) with respect to b. Hence 

K(a,P±(a))=0. (4.5) 

Thus, setting b = /3±(a) removes F(a,b) from equation (4.1). This is the key idea behind the 
"kernel method" which has been used to solve equations of this type (see [2] for example). 

Unfortunately in this case, removing the kernel reduces the recurrence to an equation con- 
taining terms F(a,ya) and F(y/3±(a), /3±(a)), which we cannot use immediately to solve for 
F(a, b). Similar situations have been studied before using the "obstinate kernel method" ([5, 6] 
for example). 

The method we use appears to be similar to the obstinate kernel method, except that instead 
of finding a finite number of pairs of values of a and b to set the kernel to zero we must use 
an infinite sequence of pairs. In this way, our "iterated kernel method" is related both to the 
kernel method and perhaps also to the iterative scheme used in [3]. 

The roots j3±(a) can be determined explicitly: 



a 



l+y 2 ± v/(l-y 2 )(l-4x|/a 2 -|/ 2 ) 



P±(a) = - ^ v ^ » " * *-L . (4.6) 

2 \ y + xa z — xy z a z 
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Define the two roots: 



Pi{a) =/3.{a) 
P_i(a)=P + (a) 



ya + 0(xy 2 a 3 ), 
a/y + 0(xy~ 2 a) 



(4.7) 
(4.8) 



as power series in a. Later we require our solution to be a formal power series in t (after setting 
x = y = t) and one can confirm that Pi(a) defines a formal power series in t. 

Since a is a variable, we are able to substitute something else for it; substituting a i— > Pi(a) 
into equation (4.5) gives 

K(p 1 (a),p 1 (p 1 (a)) = 0. (4.9) 

Hence the pair (a, b) = (/?i(a), /?i(/3i(a))) also sets the kernel to zero. We can continue in this 
way. Hence we need to define the repeated composition of Pi(a) with itself: 



(3 n (a)=f3[ n) (a)=f3 1 of3 1 o...o{3 1 (a). 



(4.10) 



Note that 



(3-i o (3i(a) = Pi o p_i(a) = a, 
P n (a) = ay n + 0(xy n+1 a 3 ). 



and 



(4.11) 
(4.12) 



There is no finite value of n such that P n = p . If we further define Po(a) = a and P- n (a) by 
composition of /3_i(a), then the functions {P n |n G Z} form an infinite group with identity p 
and inverses p n o p_ n = p . 

These observations are enough to iterate the functional equation to find a solution. Set b = 
Pi(a) in equation (4.1), and set a = P n (a) for any finite n > 0. Then since K(P n (a), P n+ i(a)) = 
0, we have 



X(p n (a),P n+ i(a))+Y(p n (a),P n+ i(a))F(p n (a),yP n (a)) 

+ Z(p n (a),P n+ i(a)) F(yP n+ i(a),p n+ i(a)) 

We can then solve this equation for F(p n (a),yP n (a)): 

~X(p n (a),p n+ i(a)) 



0. 



(4.13) 



F(P n (a),yP n (a)) 



Y(p n (a),p n+ i(a)) 
Z(p n (a),p n+ i(a)) 
Y(p n (a),p n+ i(a))_ 



F(yp n+ i(a),p n+ i(a)) 



(4.14) 



We can simplify the above by defining 

T n (a) 
X n (a) 



F{p n {a),yP n {a)) = F{yP n {a), p n {a)), 
X{P n (a),P n+ i{a)Y 



Z,„ (a) 



Y(p n (a),p n+ i(a))\ 
Z(p n (a),p n+ i(a)) 
Y(p n (a),p n+ i(a)) 



(4.15) 
(4.16) 



where we have made use of the symmetry F(a,b) = F(b,a). While this symmetry is not 
essential, it does make the solution substantially simpler. Instead of exploiting this symmetry 
we could iterate again to find F(P n (a),yP n (a)) in terms of F(P n+2 (a),yP n+ 2(a)). Indeed this is 
what is required to solve walks in the asymmetric wedge Wi (see Section 5 below). 
Equation (4.14) may be written as 



T n {a) = X n (a) + Z n (a) J" n+ i(a). 



(4.17) 
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Starting at n — 0, this can be iterated to get a series solution for To(a): 

oo n— 1 

F(a, ya) = JF (a) = ^ AT n (a) J] Z fc (a), (4.18) 

n=0 fc=0 

where we have assumed that the above sum converges (we will show that this is the case). This 
also gives F(yb, b): 

oo n— 1 

F(y6,6) = F(6,j/6) = :F„(&) = ^* n (6) J] Z,(6). (4.19) 

n=0 fc=0 

This allows us to write down the solution for F(a, b): 

'•<«•»> - ffef + B|wn««) + KJ> (i) §* (6) - «■»> 

Of course, the above "solution" still contains many complicated algebraic functions in the 
form of the (3 n (a). It is quite surprising (at least to the authors!) is that these functions can 
be drastically simplified. 



4.2 An explicit expression for /i(l,l) 

We have outlined above the iterated kernel method that we shall use to write down the gener- 
ating function fi(a, b) = F(a, b). We are primarily interested in the number of paths (and not 
the location of their endpoints), so we will actually focus on the function F(l, 1). 

We start by considering the (3 n {a) functions. It is quite surprising that while /3 n (a) is (upon 
superficial inspection for small n) very complicated, its reciprocal appears relatively simple. 
Examining equations (4.4a) and (4.6) one obtains 

1 + 1 = I±l!l. (4.21) 



f3i(a) /3_i(a) y a' 

Substituting a = /3 n _i(a) in the above, and using the group properties of (3 n leads to the 
following three term recurrence for (3 n : 

1 1+ « 2 1 1 (4.22) 



n-2 



Since (3q is the identity, and f3\ is given explicitly by ft-(a) in equation (4.6), the recurrence 
above can be iterated to get a solution for ft n (a): 

1 y(l-y 2n ) 1 y 2 (l-y 2n ~ 2 )l 



(3 n (a) y^l-y 2 )^) y n (l - y 2 ) a 

By using the expressions for X(a, b), Y(a, b) and Z(a, b) in equation (4.4) to determine X(a, b) 
and Z(a,b), one obtains 



F(a, ya) = \ Pn+1 ~f n ] ff ( ^ ~ f" ) • (4-24) 



n=0 



\Pk- yftk+i) 
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Substituting the expression for /3 n (a) given in equation (4.23) and simplifying gives: 



n+1 - yp n 
xyap n p n+1 

Pk+i ~ V(3k 

(3k - yPk+i 



,2k+l 



a Pi 
1 



and 



xya z 



xaf3\ 



Using these, one can get an explicit expression for the generating function F(a,ya): 



F(a,ya) 



xya z 



xa(3\ 



£(-i)V (n+1) 



n=0 



xya* 



xafii 



By defining 



Q(a;x,y) 



V 



xa 2 xa(3i 

the above expression for F(a, ya) can be further simplified to 

Q(a;x,yy 



F(a,ya) 



V 



J2(-i) n y n2 Q(a;x,y) n - 



(4.25) 
(4.26) 

(4.27) 
(4.28) 

(4.29) 



J n=0 



Using the a <-> b symmetry of F(a, b), we can get a similar expression for F(yb, b), and so finally 
F(a,b). 



f(a,b) 



X(a,b) Y(a,b) 
K(a,b) + K(a,b) 



n>0 



+ 



Z(a, b) 
K(a, b) 



1 



Q(b) 



V 



J2(-l) n Q(b) n y n2 (4.30) 



n>0 



We can reduce the above equation by considering only the number of walks of length n (by 
setting a = b = l,x = y = t): 

Proposition 4.1. The generating function of partially directed walks ending in a horizontal 
step in the wedge V\ is 

1-t 1-t 2 - y/(l-t 2 )(l-5t 2 ) ' 



fl ^ 1,1 ^ \-1t-t 2 \-2t-t 2 
where t counts the number of edges and 



^(-1)^3(1; M) r 



n=0 



Q(l;t,t) = (1 - 3t 2 - v/(l-t 2 )(l-5t 2 ))/2t. 
The generating function of all paths in Vi is then found using equation (3.1): 

l+t 1-t 2 - y/{l-t 2 )(l-5t 2 ) 



giO-,1) 



\-2t-t 2 



t{l - It- t 2 ) 



(4.31) 



(4.32) 



(4.33) 



n=0 



Firstly we note that F(a, ya) counts all partially directed paths in the wedge Vi whose 
last vertex ends in the line Y = —pX. Additionally we note that the generating function 
Q(a; x,y)/y counts the number of partially directed paths starting at the origin, lying on or 
above the line Y = —pX and whose last vertex lies in the line Y = —pX. Hence Q(a)/y counts 
a very similar set of paths to F(a, ya), except that the paths counted by Q are not confined by 
the line Y = pX. 

In light of the above interpretation of the function Q, we expended considerable effort to 
uncover a more direct combinatorial derivation of the alternating sum in equation (4.29). There 
appears to be some inclusion-exclusion process underlying this, but unfortunately we have not 
made progress in this respect. 
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4.3 Asymptotics for p = 1 



Where the constants are 



The asymptotics of the number of partially directed paths in the symmetric wedge with p — 1 
can be analysed by examining the singularities of the generating function g±(l, 1) in equation 
(4.33). Singularities arise either as zeros of the factor (1 — 2t — t 2 ) in equation (4.33), or as 
singularities in a/(1 — t 2 )(l — 5t 2 ), or as singularities in the series ^^L (— l) n t n Q(l;t,t) n . 

An examination of gi(l, 1) shows that it has simple poles at the solution of (1 — 2t — t 2 ) = 
0, or when t = — 1 ± \p2. We note that ^/(l — £ 2 )(1 — 5t 2 ) has branch-points (square root 
singularities) at t = ±1 and again at t = ±l/y/E. The series ^2^ =0 (— l) n t" 2 Q(l; t, t) n is a 
Jacobi ^-function and it is convergent inside the unit circle except at singularities of Q(l;t,t); 
that is, when t = ±1/ v5. 

The dominant singularity is the simple pole at y/2 — 1 , while the next sub-dominant con- 
tributions to the asymptotics will be given by the singularities at t = ±1/\^E. These two 
sub-dominant singularities will give a parity effect. The contributions from these singularities 
allow us to write down the asymptotic form of t> n .i- 

Proposition 4.2. The number of paths in the wedge V± is asymptotic to 

<V = A)(l + V2) + (A, + (-1) M 2 + 0(l/n)j . (4.34) 

A = 0.277309853486031 18827 . . . , (4.35a) 
A l = 3.71410486533662324953 . . . , (4.35b) 
A 2 = 0.20697997020804157910 .... (4.35c) 

We note that the constants were derived by expanding the expression for ^(1,1) about 
t = \f2 — 1 and t = ±1/^/5 (or rather the first 40 or so terms of the sum). These were then 
checked using both Bruno Salvy's gdev package for Maple [24] and by direct examination of 
v n> \ for n < 1000. The above formula is quite precise and it correctly estimates fio,i, "^20,1, ^30,1 
and ^40,1 to within 7%, 1%, 0.2% and 0.06% respectively. 

Note that the above result implies that walks in the wedge V p have the same dominant 
asymptotic behaviour as walks with no bounding wedge (see equation (1.4)). Since the number 
of walks in any wedge V p for 1 < p < 00 is bounded between the number of walks in Vi and 
partially directed walks with no bounding wedge, we have the following result: 

Corollary 4.3. The number of partially directed walks in the wedge V p , cjfi obeys the following 
inequality 

(p) 

0.2773 . . . < lim — _ < (1 + V2)/2 = 1.2071 . . . (4.36) 

n-+oo (1 _|_ y/2) n 

for any 1 < p < 00. 



5 Partially Directed Paths in the Asymmetric Wedge 

In this section we turn our attention to the model in Figure 2 (right). The partially directed 
path is confined to an asymmetric wedge, W p , and its generating function does not have the 
a <-» b symmetry we have exploited in solving for fi(a, b) in the previous section. 

We proceed by examining the generating function of walks that end in a horizontal step. 
The functional equation for these walks is given in Proposition 3.2 and we can arrange equation 
(3.5) in kernel form: 

K(a, b) h p (a, b) = X(a, b) + Y(a, b) h p (a, ya) + Z(a, b) h p (yb, b) (5.1) 
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where 



K(a, b) = (b- ya){a - yb)(l - xa p ) - xya p {a 2 + b 2 - 2yab), (5.2a) 

X{a,b) = (b-ya){a-yb), (5.2b) 

Y(a, b) = -xya p+1 (a - yb), (5.2c) 

Z(a, b) = -xya p b(b - ya). (5. 2d) 

This functional equation is very similar to that of the symmetric wedge given in equation (4.1). 
However we no longer have a b symmetry and this means that we have to work quite a bit 
harder and we concentrate only on the case p = 1. We will write hi(a,b) = H(a,b) for the 
remainder of this section. 



5.1 Solving for H(a, b) when p = 1 

For the remainder of this section we concentrate on the case p = 1 and walks in the 45° wedge 
Wi. The kernel K(a, b) (given in equation (5.2)) is no longer symmetric in a and b, nor is the 
desired generating function H(a, b). In order to repeat the iterated kernel method as described 
in Section 4.1 we must now consider the solutions of the kernel as functions of a and b. These 
solutions are defined by K(a,/3(a)) = and K(a(b),b) = 0: 



At (a) 



a 

2y L 



and 



y 2 - x(l - y 2 )a ± - y 2 )((l - xa) 2 - y 2 {l + xa) 
l+y 2 ±^/(l-y 2 )(l-y 2 -Axyb) 



a±(b) = 2 



y + x(l — y 2 )b 



(5.3) 



(5.4) 



One can confirm that a_(6) and /3_(a) both define formal power series in b and a (respectively). 
Additionally these same choices (when x = y = t) also define formal power series in t — which 
we will require for our solution. Write these as a\(b) and (3i(b), and the other roots as a-i(b) 
and (3-\(a). 

In Section 4.1 we considered composing j3{a) with itself, however due to the asymmetry of 
the kernel we now need to consider mixed compositions /3(a(b)) and a ((3(a)). Indeed, we find 
that 



a±iC&Fi(a)) = a, 
f3 ±1 (a Tl (b))=b. 



(5.5a) 
(5.5b) 



We will need the function 7(a) = «i(/3i(a)), and define its nested composition by 7„(a) = 
7(7„_i(a)) with 70(a) = a. Note that 

ln (a)=y 2n a + 0(xy 2n a 2 ). (5.6) 

We can now repeat the iterated kernel method in the new asymmetric setting. Setting 
b = Pi(a) in equation (5.1) gives 

= X{a, /3i(a)) +Y(a, i3 1 {a))H(a,ya) + Z(a, /3i(a))if(y/3i(a), Pi{a)). (5.7) 

Since there is apparently not a simple relation between H(by,b) and H(b,by), this equation 
cannot be iterated to find a solution. Instead, it turns out that the other roots of the kernel 
must be considered as well. 
Setting a = «i(6) gives: 



= X{ ai (b), b) + Y( ai (b), b)H( ai (b), y ai (b)) + Z{ ai {b), b)H(yb, b). 



(5.8) 
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Now set b = Pi (a) in the above equation 



= X{ 1 (a),(3 1 (a))+Y( ll (a),p 1 (a))H( ll {a),y ll (a)) 

+ Z( ll (a),p 1 (a))H( y p 1 (a),p 1 (a)). (5.9) 

We can now eliminate H(yPi(a), Pi(a)) between equations (5.7) and (5.9) and solve for if (a, ya). 
This gives 



H(a,ya) = - 



XjaMa)} 



+ 



Z(a,p 1 (a)) 
Y{a^x{a)) 



X( 7l (a),/3 1 (a)) 
Z( 7 i(a), &(<*)) 



+ 



Z(a,p 1 (a)) 



y(7i(a),A(a)) 
Z( 7l (a),ft(a)) 



#( 7l (a),?/ 7l (a)). (5.10) 



We can now iterate the above equation by substituting a = 7n _i(a). Define 

X( 7n (a),/5i( 7n (a))) P{ln)-yin 



^n(a) 
ZJa) 



Y( ln (a),p 1 ( ln (a))) 
Z(7 n (a),/3i(7 n (a))) 
F( 7n (a),/3i( 7n (a))) 
-X"(7n+i(a),^i(7n(a))) 

Z(7n+l0),/?l(7n(a))) 



/?(7n) / /?( 7 n) - Vln 



In \7n - 
7n+ l - ?//?( 7n ) 
£?/ 7 n+l/?( 7 n) ' 



A ( 1 = F (^+ 1 ( Q )^ 1 ^( a ))) 7n+l / 7r»+l ~ l/^(7n) \ 

Z( 7n+1 (a)), A(7n(a))) /3(7n) W) - Vln+J ' 



C„(a) = y n (a)An(a) 



And further define 

B„(a) = X n {a) + 3> n (a)Z n (a), 
Equation (5.10) now becomes: 

#( 7n (a), j/7„(a)) = B„ + C n F( 7n+ i(a), y 7n+1 (a)). 
We obtain a solution for H(a, ya) by iterating the above equation 



n-l 



H(a,ya) = B + C B 1 +C C 1 B 2 



(5.11a) 
(5.11b) 
(5.11c) 
(5. lid) 

(5.12) 
(5.13) 

(5.14) 



n=0 



m=0 



As was the case for the symmetric wedge, we are able to simplify the above expression by 
rewriting 7n (a) in terms of the original kernel roots and thereby rewrite the expressions for B n 
and C n . 

In order to find H(a,b) from equation (5.1), we need both H(a,ya) andH(yb,b). Equa- 
tion (5.8) gives H(b,yb) in terms of H(a\(b),yai(b)): 



H{yb,b) 



X( ai (b),b) Y( ai (b),b) 



Z( ai (b),b) Z( ai (b),b) 



H(a 1 (b),ya 1 (b)) 



(5.15) 



So using the above expressions for H(a, ya) and H(yb, b) we have, at least in principle, a solution 
for H(a,b): 



K(a,b) Z(a 1 (b),b)K(a,b) K(a,b)' 



Z(a,b)Y(a 1 (b),b) 
Z(an(b),b)K(a,b) 



^(ai(6),yai(6)) (5.16) 
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Of course, we would like to be able to simplify the above expression. In particular we would like 
to rewrite 7 n (a) and 7 n (ai(6)) in terms of simpler functions, as we did for (3 n (a) in Section 4.2. 

Interestingly enough, it is possible to determine H(a,ya) in equation (5.14) by inspection 
of the terms in this expression. Putting x = y = t, one obtains 



if (a, ta) 



£ 

n=0 



£2(n+l) 2 -3 



a-pt- apt 2 + apt 2n+2 



_a(l + p)t 2n ■ - (3(a + t 2 ^ 1 ) 
a-pt- apt 2 - pt 4n+1 + a(l + p)t 4n+2 



n 

m=0 



+ (a(l - /3)t 2 + a(l + /?)t 2 « +2 ) 

a(l + /3)t 2m - ff(a + t 2 '"- 1 ) 
a-pt- apt 2 + a/% 2m + 2 



(5.17) 



From this expression one may determine H(bt,t) from equation (5.8), and thus an expression 
for H(a,b). While the resulting expression gives a series expansion for the numbers of paths, 
it is not very useful because it is so complex. In the next section we proceed by simplifying 
expressions for the compositions of the a's and /3's above; ultimately this will lead to a simpler 
expression for H(a,b). 



5.2 Simplifying things 

In much the same way as for the symmetric case, we can find simple expressions for the l/7 n (a) 
in terms of the original kernel roots. Consideration of the kernel and its roots gives: 

1 1 1 + y 2 . . 

(5.18a) 



ct_i(6) a + i(b) yb 

1 ■ 1 1 + ^ (5.18b) 



P-i(a) P+i(a) ya y 

Since certain compositions of a and P give the identity (see equations (5.5)), we have the 
additional relations: 

(5.19a) 



aii(Pi(a)) 71(a) yPi(a) a 

1 _ 1 + y 2 1 x(l -y 2 ) 
Px(a x {b)) yax{b) b y 



(5.19b) 



We note that the last term in equation (5.19b) means that the resulting expressions for 7 n (o) 
are more complicated than those for P n (a) for the symmetric case (see equation (4.23)); this in 
turn leads to a significantly more complicated solution. 

Setting a = 7„_i(a) and b = /5i(7„_i(a)) in the above two equations give: 

l l + y 2 l , 

(5.20a) 



7n(o) yPi(ln-i(a)) 7n-i(a)' 
1 l+y 2 1 x(l-y 2 



Pihn(a)) yin(a) /?i(7 n _i(a)) y 



(5.20b) 
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These equations can be solved: 

1 l-y 4n l-y 4n - 2 x{l - y 2n ){\ - y 2n ~ 2 ) 



7 n(a) y*«-\l-y 2 )Px{a) y 2n ~ 2 (l-y 2 )a y 2 ^ 2 (l-y 2 ) 

= T^—i H 1 + y 2 ) + Q( a ) y 2n + y~ 2n ) > ( 5 - 21a ) 

1 — y z v 7 

1 1 _ y 4n+2 1 _ y in ^ _ y 2r^2 



A( 7 n(a)) ^(l-y 2 )/^) _ y 2 )a y 2n-i {1 _ y 2y 

(2xy + yQ(a) y 2n + Q(a) y~ 2n ) . (5.21b) 



i - y 



2 



where we have used 



Q(a) = --TrU-x Q{a) = -^—- y --xy (5.22) 
a pi(a) Pi (a) a 

Note that Q(o)Q(o) = x 2 y. In fact we can reduce the above expressions for 7„ and (3{ r j n ) even 
further using this fact: 

1 _ (x + y 2 »- 2 Q)(z + y 2 »Q) ^ 



7 „(a) </ 2 - 2 (l-y 2 )Q 
1 (x + y 2n Q) 2 



^(7n(a)) ^-Hl-y 2 )^ 



(5.23b) 



The above then lead to the following expressions that will be useful in writing down our 
solution: 

1 V ... , _.2n 



= (x + ^Q), (5.24a) 
7n(a) Pi(7«(a)) 

mM)-^^^' (5 - 24b) 

ASw)-dw^ (l+!/2 " Q) ' (5 - 24c) 

1 y x 



(x + y 2n Q) (5.24d) 



7„+i(a) A(7n(a)) y 2n( 5 

This in turn lets us write 

g(7n) ~ 2/7n _ Z/ 2 "" 1 

7n - y/3(ln) X 
7 n+l ~ y£(jn) = y 2 " +1 
/%n) - y7n+l « 

(5.25c) 

where we have made use of the fact that Q = x 2 y/Q. Hence C n = y n A n can now be written as 

C = ^±1 ( P^hl ~ Vln \ ( 7n+l -yf3(ln) \ = 7n+l _ 2/___g2 ^ 2 g) 



Q (5.25a) 
Q (5.25b) 



7n \ 7fi 

In a similar way we find that 

Piln) ~ Vln = y(x + y 2n ~ 2 Q) (5.27a) 



7n 

y 2 (x + y 2n ~ 2 Q) (5.27b) 



7n+l - y/%n) 



7n7n+l 
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which allows us to also simplify the expression for B n — X n + y n Z n : 
xy \ % J x 

y nZri = _L h^-v^ )( ^-^A = _£*L (x+1 ?~* Q) (5 . 28b) 

Xy \ 7n7n+l J\ln-yp{ln)J X 2 

B n = X n + y n Zn = ^{x + y 2n ~ 2 Q) {x - y 2n Q) (5.28c) 

Substituting the above into the expression for H(a,ya) in equation (5.14) gives: 

oo n— 1 

H(a,ya) = ^B n {a) JJ C m {a). 

n=0 m=0 

= E A (* + y 2n ~ 2 Q) ( x - y 2n Q) — (-) ^ y 2n[n ~ 1] 

X CL \ X / 

•n— n \ ' 



n=0 



which is a significant simplification of equation (5.17). We can now substitute this into equa- 
tion (5.16) to obtain H(a,b). This requires us to compute H(ai(b), ya\{b)) from the above 
expression. Let 

P(b) = Q(ai(6)) = ^ (l - 2xy6 - y 2 + ^(1 - y 2 )(l - 4^6 - y 2 )) (5.30) 

then we have 

tt( \ (l-y 2 )P ^ (x-y 2n P) (P\ 2n 2n > 

F( ^ a) = ^^^P)UJ y (5 ' 31) 

Below we give the length generating function (when x = t,y = t,a = l,b = 1). 

Proposition 5.1. The generating function of partially directed walks ending in a horizontal 
step in the wedge Wi is 



i) = — ™ — 7^ 



2(l-2t-t 2 ) 

c (i - 1 2 ) f, (i - t 2 *-^) /qn 2n 

^t 2 (l -It- t 2 ) ^ (1 + t^-iQ) V * 

v ' n=0 



| a-^) f (l-^P) ^ ^ 



1 - 2t-t 2 ^ (1 +t 2n ~ 1 P) V t 
n=0 v 7 v 

where t counts the number of edges and 

Q(l;t,t) = (1 -t-t 2 -t 3 - y/(l-t*)(l - 2t-t 2 ))/2 (5.33a) 
P(l; t, t) = (1 - 3t 2 - V(l - t 2 )(l - 5t 2 ))/2t (5.33b) 

The generating function of all paths inWi is f/ien fci(l, 1; t, t) = (hi(l,l;t,t) — l)/t. 

As was the case for the symmetric wedge, the functions P and Q that make up our expression 
for H have combinatorial interpretations in terms of partially directed walks bounded by a single 
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line. Let B_(x,y) be the generating function of walks that end with a horizontal step, start 
and end on the line Y = and stay on or above that same line. Then 

Q(a;x,y)=xy 2 (B4ax,y)-l). (5.34) 

Similarly let B/(x, y) be the generating function of walks that end with a horizontal step, start 
and end on the line Y = X and stay on or above that same line. Then 

P(b;x,y)=xy 2 (B f (bx,y)-l). (5.35) 

Again we would like to find a more direct combinatorial derivation of the generating functions 
H(l, 1) and H(l,t). We have been unable to do so. 

5.3 Asymptotics for p = 1 

Before we study the asymptotics of walks in the wedge Wi, let us compute the number of walks 
lying on or above the line y = 0. 

Lemma 5.2. The generating function of partially directed walks lying on or above the line 
y = is 



-1 + z + 3z 2 + z 3 - - z 4 )(l -2z- z 2 ] 
2z 2 (z 2 -2z-l) 
The number of these walks is asymptotic to 



(5.36) 



7 + 5^(l + v^)" (1 + 0(1/n)) 



2tx ^/n 

Hence the number of walks in the wedge Wi must also be 0((\/2 + l) n /\fn). 

Proof. One can derive a functional equation for the generating function of such walks which 
can be solved using the kernel method. The asymptotics can then be computed by analysing 
the dominant singularity at t = \/2 — 1. The last result follows since the number of walks in 
the wedge Wi cannot exceed the number of walks lying above y = 0. □ 
As was the case for walks in the symmetric wedge, we analyse the asymptotics of partially 
directed walks in the asymmetric wedge Wi by singularity analysis. Let us split the expression 
given in equation (5.32) into 3 pieces and study their dominant singularities: 



to - (l-*)W(l-* 2 )(l-5t 2 ) f5 3ga) 

Pl ~ 2(l-2t-f) { } 



-t 2 ) ~ (Q^ 

n=0 

[l-t 2 ) ^ (l-t 2n ~ l p) (p 



P3 



1 - 2t - t 2 ^ (1 + t 2n - x P) V t 

n=0 K 



T ) L ~\ " - t 2 - 2 (5.38c) 



Let us treat the asymptotics of each of these functions separately. 
Lemma 5.3. The coefficients ofp±(t) are asymptotic to 



[t n ]Pi 



- ^ • ((2 + V5) - (-1)"(VS - 2)) • ■ (1 + 0(n~ 1 )) (5.39) 
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Proof. The generating function p\ appears to have 6 singularities: 2 simple poles from the zeros 
of the denominators and four square-root singularities at t — ±1, ±1/ Closer analysis shows 
that there are no singularities at the zeros of the denominator and that generating function is 
dominated by the singularities at t = ±1/\^E. Analysis (by the techniques in [14] ) of these 
singularities leads to the above expression. □ 
Before we can study the asymptotics of P2(t) and ps(t) we need the following lemma about 
the location of the zeros of 1 + Qt k and 1 + Pt h . 

Lemma 5.4. Fork = —1,0, 1,2, . . ., the functions 1 + Q(t)t k and \ + P{t)t k are not zero within 
the disk \t\ < 1/2. 

Proof. The function Q(t) satisfies Q 2 — (1 — t — t 2 — t 3 )Q + 1 4 = 0, and so h = Qt k satisfies: 

h 2 - (1 - t - t 2 - t 3 )ht k + t 2k+A . (5.40) 

Now if 1 + Qt k = we have h = — 1 and so 

l + (l-t-t 2 -t 3 )t k + t 2k+4 = 0. (5.41) 

For \t\ < 1/2 we can bound |1 - t - t 2 - t 3 \ < (1 + \ + \ + |) < 2, and therefore 

|(1 - t - t 2 - t 3 )t k + t 2k+4 \ < \2t k + t 2k+i \ < 3\t\ k <3-2~ k (5.42) 

Hence for k > 2, the above quantity is less than 1 and so equation (5.41) cannot be satisfied. 
It remains to check the cases k = —1, 0, 1. In these cases we can solve equation (5.41) directly 
and verify that t lies outside |t| < 1/2. 

The argument for P(t) follows a similar line. The function h = Pt k satisfies: 

h 2 -t k ~\l-3t 2 )h + t 2k+2 = 0. (5.43) 

Hence if h — — 1 we have 

1 + (1 - 3t 2 )t k - 1 + t 2k+2 = (5.44) 
For \t\ < 1/2 we can bound |1 — 3t 2 | < 2 and so 

|(1 - 3t 2 )t k ^ + t 2k+2 \ < 3\t\ k ~ l (5.45) 

Hence for k > 3, equation (5.44) cannot be satisfied for \t\ < 1/2. For k = —1,0, 1,2, we can 
check equation (5.44) directly and verify that the zeros do not lie inside |t| < 1/2. 

Note that when k = 0, equation (5.44) has a solution t = \[2 — 1, however this point 
corresponds to the other branch of P being +1. 

□ 

We can now move onto the asymptotics of P2{t) and P3(t). 
Lemma 5.5. The coefficients of p^it) are asymptotic to 

[ mKJ V2 1 + (V2 - l) 2k+lK ' VV / 

= (0.31096381899209832 . . .) (V2 + l) n + O ( (V^J ") (5.46) 

Proof. The function p 3 has a simple pole from the zero of the denominator of prefactor. There 
are also square-root singularities in P{t) at t — ±1, ±l/\/5, simple poles when 1 + Pt 2n ~ l = 
and a natural boundary at \t\ = 1 coming from the ^-function like structure of the sum. Of 
these singularities, the simple pole dominates, followed by the singularities at ±l/-\/5- The 
simple pole and its residue give the dominant asymptotics and the square-root singularities 
give the 0(5 n / 2 ) corrections. □ 
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Lemma 5.6. The coefficients of p2{t) are asymptotic to 



[t> 2 (t) = - (v ^+ 1)w y 1 ~ l)2k+ \ v~2 - i) 2fc2+2fc (i + o(i)) 

= -(0.31096381899209832 ...)(V2 + l) n (1 + o(l)) (5.47) 

The dominant term is equal in magnitude, but opposite in sign, to the dominant term in the 
asymptotics ofps(t). 

Proof. The singularities of the function P2{t) arise from singularities of the prefactor, the singu- 
larities of Q(t), the zeros of 1 + Qt 2n ~ 1 and the natural boundary at |t| = 1 from the ^-function 
structure of the sum. Of these, the simple pole of the prefactor at t = \/2 — 1 and the singu- 
larities arising from Q(t) at the same point, dominate the asymptotics. 

The contribution from the simple pole may be computed by finding its residue. We note 
that P( v / 2 - 1) = Q(y/2 - 1) = 3 - 2y/2, and this means that the residue is in fact equal in 
magnitude, but opposite in sign, to that computed for ps(t). 

□ 

We note that if the dominant asymptotics of p2(t) and ps(t) must cancel each other. Oth- 
erwise, the number of walks in this wedge is ~ (\/2 — l) n which would contradict Lemma 5.2. 

Lemma 5.7. The k th summand of P2{t) is 

v (t)-( c (1 " t2) ^- t2k ^) (Q\ 2 \A (5 48) 

P2 > k{t) -{- Q ti(l-2t-t*)(l+ t^Q) [l ) 1 ) (5 - 48) 
so that p2(t) = J2Pz,k(f)- The coefficient off 1 in p 2 ±(t) is asymptotic to 



r w( t)= -( 1+ ^r.^( ;;^ : ;p (i + v2)- 



J — 

Tin 




{2k + 1)(1 - (V2 - l) 4fc+2 ) - (V2 - l) 2fc+1 



;i + ( v / 2-i) 2fc+1 ) 2 
1 



+0 



(l + V2)- 2k2 - 2k - 5 / 2 (5.49) 



Proof. The result follows from standard singularity analysis ([14]) of P2,k{t)- D 

Remark. It is unfortunately the case that we have been unable to proceed completely rigor- 
ously from this point. In particular, we have been unable to obtain uniform bounds on the error 
terms in the above asymptotic expressions. On the basis of numerical testing, we think that 
the that the error term is 0(k 3 / \/n). If this is the case, then one can sum the contributions of 
the individual P2,k{t) to obtain the asymptotics of coefficients of P2{t). 

We believe that the expressions that follow are indeed exact, if not completely rigorous. 

Assuming that the asymptotic expression in the previous lemma has a uniform error bound, 
so that we may sum the contribution to the individual P2 k {t), we find that 



[t n }p 2 (t) = (1 + V2) n [ - 0.31096381899209832 . . . 

0.090584741026764287. 



0(1/Vn3)) (5.50) 
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where the constant 0.31 ... is the constant that appears in both Lemmas 5.5 and 5.6. So adding 
together the contributions from the pi we obtain 



MMi, i) = (i + V2)" ( 090584741 J 6764287 ^ + o(i/V5i)) 



(5.51) 



The generating function hi(l,l) enumerates walks that end in a horizontal step, and so we 
obtain the number of walks in Wi ending with any step by multiplying this expression by a 
factor of (1 + y/[2)) ( since the generating functions differ by a factor of t): 

rM l, !) = (! + V2)» ^218693916^4303177... + Q(y ^ ^ 

We have confirmed this numerically using the first 1000 terms in the series expansion of fci(l, 1). 
Additionally we have check the above expression using Bruno Salvy's gdev package for Maple 
[24]. 

Remark. We note that if the above result is indeed made rigorous then we have a result 
analogous to Corollary 4.3. The number of walks in any wedge W p for 1 < p < oo is bounded 
between the number of walks in Wi and partially directed walks inside the first quadrant (see 
Lemma 5.2. Hence the number of partially directed walks in the wedge W p , cif 1 obeys the 
following inequality 

c>y/ 2 




0.21869 . . . < lim — - — — < \ — — — = 1.496489 . . . (5.53) 
- ™ (i + y/2)r 



for any 1 < p < oo. 



6 Conclusions 

In this paper we have proved that partially directed paths in the wedges V p and W p all grow 
with the same exponential growth rate 1 + independent of p. Additionally we have found 
generating functions for partially directed paths in the symmetric wedge Vi and the asymmetric 
wedge Wi, using a variation of the kernel method. From these generating functions we have 
computed the asymptotics of the number of paths in both of these wedges. 

Curiously the number of paths in the symmetric wedge, V\, has the same leading asymptotic 
behaviour as partially directed paths with no bounding wedge. Similarly the number of paths 
in the asymmetric wedge, Wi, has the same leading asymptotic behaviour as partially directed 
paths above the line Y — 0. Because of this, we are able to determine the leading asymptotic 
behaviour of paths in the wedges V p and W p for all p > 1. 
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